Two-point density correlations of quasicondensates in free expansion 
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We measure the two-point density correlation function of freely expanding quasicondensates in the 
weakly interacting quasi-one-dimensional (ID) regime. While initially suppressed in the trap, den- 
sity fluctuations emerge gradually during expansion as a result of initial phase fluctuations present 
in the trapped quasicondensate. Asymptotically, they are governed by the thermal coherence length 
of the system. Our measurements take place in an intermediate regime where density correlations 
are related to near-field diffraction effects and anomalous correlations play an important role. Com- 
parison with a recent theoretical approach described by Imambekov et al. yields good agreement 
with our experimental results and shows that density correlations can be used for thermometry of 
quasicondensates . 

PACS numbers: 03.75.Hh, 63.22.-m, 71.45.Gm 
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Interacting Bose gases confined in strongly elongated 
traps exhibit unique properties and quantum phases re- 
lated to the one-dimensional (ID) character of the under- 
lying physics [1-3]. They also represent one of the few 
complex many-body systems which allow a direct com- 
parison with exact and often analytical theoretical de- 
scription [4j. For such gases, the smooth transition from a 
fully decoherent system to a finite-size Bose-Einstein con- 
densate is characterized by an intermediate quasiconden- 
sate regime where density fluctuations are suppressed [I[ 
whereas axial phonon-likc excitations are thermally pop- 
ulated. As a result the gas displays ID phase fluctuations 
along its axial direction, strongly affecting its coherence 
properties. 

Early experimental evidences for this regime have been 
obtained in expansion experiments where phase fluctua- 
tions of the trapped quasicondensate turn into density 
fluctuations @. Momentum Bragg spectros copy Q and 
more recently matter wave interferometry [a-llO] have al- 
lowed one to investigate in depth the first-order coherence 
properties of such systems. More involved interferome- 
try schemes have been used to measure the second order 
correlation function of trapped quasicondensates 
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In this letter we report on the experimental measure- 
ment of the two-point correlation function of quasicon- 
densates released from a tight atom-chip wave guide. The 
theoretical description of quasicondensate expansion is 
motivated by the possibility of a unique mapping of the 
emerging density correlations back to the trapped sys- 
tem [13j. In previous experiments, this was complicated 
by interaction-induced hydrodynamic effects during ex- 
pansion [rj.ll3|. A complete description of density correla- 
tions then requires involved numerical calculations. Our 
samples are sufficiently deep in the ID regime so that the 



radial dynamics can be neglected and the expansion can 
be considered as collisionless. This allows us to apply 
the theory of Imambekov et al. [13j and directly relate 
the density correlations in expansion to the correlation 
length of phase fluctuations present in the trapped sys- 
tem. Our measurements take place in a near-field regime 
where the shape of the two-point correlation function re- 
flects the transition from the trapped system (suppressed 
density fluctuations) to the asymptotic case (exponential 
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FIG. 1: (a) Linear density profile for a short expansion time 
(2 ms), where density fluctuations have not yet developed, and 
fitted density profile [l|. (b) Linear density profile after 10 
ms of free expansion, showing high contrast density fluctua- 
tions, where the averaged density profile (gray dashed line) is 
smooth, (c) Example images of density fluctuations emerging 
during free expansion. 
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decay on a length scale given by the in-trap coherence 
length). Extending on the work with elongated three- 
dimensional (3D) Bose-Einstein condensates [l4| we use 
the amplitude of density fluctuations of quasicondensates 
for thermometry in a temperature regime where conven- 
tional methods based on the visibility of thermal atoms 
fail. 

The starting point for our experiments is a thermal 
cloud of 87 Rb atoms in the \F = l,mp = —1) state pre- 
pared in a strongly anisotropic Ioffc-Pritchard-type wire 
trap on a multilayer atom-chip [15']. The measured trap 
frequencies are io r — 2tt x 2.3 kHz in the radial and 
uj a = 2tt x 19 Hz in the axial direction. We use forced RF 
evaporation to cool the gas into the quantum degenerate 
regime, the final temperature is adjusted by limiting the 
effective trap depth to (1 — 10) huj r . To ensure thermal 
equilibrium, this situation is maintained for several hun- 
dreds of milliseconds. We adjust an additional hold-time 
to keep an atom number of rs 3500 in all presented exper- 
iments, making use of technical losses. This corresponds 
to a peak atomic density of ?iid(0) = 90 /im -1 and a 
chemical potential of fi/h — 2ir x 1.6kHz= 0.7 u> r [l2j ]. 
For temperatures T < huj r /kB = 110 nK we realize 
the regime of weakly interacting quasi-lD Bose gases 
([/,, IcbT < fou r ) where the radial degrees of motion have 
negligible influence on the system 20] . In this regime the 
quasicondensate is characterized by thermal phase fluc- 
tuations with a correlation length At = 2h 2 niY)/mkBT, 
m being the mass of the atoms. For our experiment At 
ranges from 0.25i? a to 0.5i? a where 2R a = 64/Ltm is the 
axial size of the gas. 

After preparation of the quasicondensate, the confin- 
ing atom-chip potential is suddenly switched off (within 
less than 10/xs) and the atomic density is imaged after 
an expansion time i exp using standard absorption imag- 



ing [21] . During the switchoff, the expansion changes 
from hydrodynamic to collisionless on a time scale on 
the order of cj~ . During this time, short-range fluctua- 
tions up to the size of 5x — £h(i/hu) r , £/> = h/^Jmp, being 
the healing length, will be smeared, and are not visible 
in the correlation function [l3j]. Since this length scale 
is on the order of 100 nm it cannot be resolved with our 
imaging system. For the visible density fluctuations we 
can therefore assume ballistic expansion. 

Typical experimental images are depicted in Fig. [1] 
Starting from a smooth in situ density profile, the atomic 
density develops density fluctuations in the course of the 
expansion, reflecting the thermal phase fluctuations ini- 
tially pre sent in the quasicondensate. In contrast to other 
work [l6l Il7j |. density fluctuations in the initial sample 
(e.g., due to corrugation of the trapping potential) can 
be neglected here. In expansion, each individual real- 
ization exhibits different density fluctuations. The mean 
image, however, shows a smooth profile. 

The theoretical model by Imambekov et al. [l3| de- 
scribes the axial shape of the density correlation func- 
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FIG. 2: Two-point correlation function of an ultracold ID 
Bose gas for different times t cxp after the release from the trap 
in the case of weak interactions (quasicondensate limit). The 
units on both axes are dimensionless. In the trap (t cxp — 0) 



g2(x) ~ 1, the asymptotic limit 
l + ffltp[-2[x|/(AT)]. 
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tion g2(x,t exp ) of a uniform ID Bose g clS clS Si function 
of expansion time, ID density and temperature (com- 
pare fig. [2]). To allow for comparison with our experi- 
mental finite-size quasicondensate, we determine an av- 
eraged two-point correlation function c/2 defined as fol- 
lows : We first compute the autocorrelation function for 
a single integrated density profile J nm(u)nm(u + x)du 
and obtain Gzix) — (j niD(u)niu('u + x)du) by aver- 
aging over many (several hundred) repetitions. The 
function gi is then obtained by normalizing with the 
autocorrelation function of the mean density profile 
/ (niT>(u)) (ni D (u + x)) du. 

Figure [3] shows experimentally obtained two-point den- 
sity correlation functions for two different temperatures 



after an expansion time t e 



10 ms. We observe corre- 



lation well in excess of 1 (bunching) for small relative dis- 
tances, as expected for an only partially coherent source 
such as a phase fluctuating (and hence multimode) quasi- 
condensate. Furthermore, we find minimum correlations 
below 1 at finite distance, which is a signature of the 
transition from the trapped gas g2(x,t cxp = 0) « 1 to 
the far-field limit 32(2^ ^oxp 00) ~ 1 + exp(— 2 |x| /At) 
throughout the free expansion. 

This behavior can be understood in terms of probabil- 
ity flow in the relative distance space towards x = re- 
sulting in the deficit of probability to find a second atom 
within a certain distance range. As t exp increases, this 
range spreads toward larger The depth of the min- 
imum decreases as i C xp — > 00, and finally the minimum 
vanishes for all finite |x|. Note that the surface below 
the g2(x, t C xp) function is maintained for all expansion 
times. Alternatively, this behavior may be described as 
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near-field diffraction in terms of a matter wave analogous 
to the Talbot effect [l~3] . The two-particle problem is re- 
duced to a single particle one by elimination of the center- 
of-mass motion. The random phase of the matter wave 
in the relative coordinate is even due to Bose-Einstein 
statistics, which leads to a g% (x) peak formation always 
at x = 0. As illustrated in Fig. [2] the oscillatory behav- 
ior will only vanish for expansion times on the order of 
seconds. 

To compare theory and data, the finite resolution of 
the imaging system has to be taken into account. Typ- 
ical structure sizes in the density fluctuations are on 
the order of the system's point spread function (PSF), 
which hinders the observation of smaller features. The 
PSF can be determined, starting from a description in 



Fourier space [18|, which includes effects of diffraction, 



pixel sampling and finite extent of the cloud along the 
optical axis. It is then approximated by a Gaussian 
curve with a half width at half maximum (HWHM) of 
a i = 3.9 /im. For two-point correlation data, this value 
becomes 02 = V^ai. The PSF is finally convolved with 
the theoretical prediction for the correlation function 
g2{x). Figure [3] compares a measured 32(2:) function to 
a theoretical calculation based on the experimental pa- 
rameters and the imaging resolution. We find excellent 
agreement without any adjustable parameters. 

As the expansion time and the atomic ID density are 
known to high precision, the density correlation function 
can be used to determine the temperature of the system 
in cases where this information cannot be retrieved other- 
wise. In the interesting temperature range from 20 nK to 
the critical temperature T c ss 350 nK, the density corre- 
lation function at zero distance 32(0, t oxp ) is a monotone 
function of T that can be calculated in a straightforward 
fashion from [l3| and can hence be used for thermometry 
(compare Fig. 2]). 

Before comparing data to theory, again technical as- 
pects of the imaging system have to be considered. 
The noise in absorption images is typically composed of 
atomic and photonic shot noise, and technical detector 
noise. We study the noise characteristic of our system by 
evaluating images integrated along the quasicondensate 
axis with and without atomic signal. It shows that the 
dominating contribution to noise is technical. Atomic 
shot noise can be neglected since the main contribution 
to the variance on a single pixel arises from the density 
fluctuations themselves. The technical noise contribution 
is uncorrelated and therefore contributes to 32(0, i e xp) 
only. To account for this noise in the analysis, we use an 
empty region of each absorption picture and integrate it 
in the same direction as the region of interest containing 
the atomic signal. Then we calculate the variance of the 
empty region and subtract this value from the autocorre- 
lation of the corresponding integrated region of interest 
before averaging and normalization. 

To evaluate whether measurements of density correla- 




FIG. 3: Black dots (open diamonds): Measured correla- 
tion function c/2(x) for 150nK (60nK) [ i.e., M exp /m\T = 
0.18 (0.03)]. The solid line (dashed line) shows the according 
theoretical prediction convolved with the imaging PSF. (In- 
set) Solid line (dashed line): corresponding theory curves not 
considering the imaging system. 



tion functions can be used for quantitative comparison, 
we need an alternative temperature calibration for our 
samples. For the quasi- ID regime that is currently ac- 
cessible in our setup, the condensate fraction is between 
70% and 90%. For higher temperatures (T > T c /2) a 
more pronounced thermal fraction is present and the tem- 
perature of the sample can be determined by fitting a 
Bose function to the tails of the density profile. For low 
temperatures, we perform a similar measurement using 
a special fluorescence imaging scheme [l8j], which allows 
us to detect small thermal fractions down to 10% in a 
single measurement 22]. To gauge the temperatures of 
our samples, we relate the temperature of a single mea- 
surement to the final value of the RF frequency used for 
evaporative cooling. We find a linear dependence over a 
range 4 T c to 0.2 T c . 

In fig. Hwe compare experimental values of 32(2^ = 



(Me 



obtained for different temperatures to the theo- 



retical prediction based on |13| and find agreement within 
the error bars. Note that we expect the theory to devi- 
ate from the data for temperatures above 110 nK as the 
system becomes increasingly three-dimensional. Temper- 
atures obtained from 32(2; = 0, t eKp ) compared to temper- 
atures derived from the calibration described previously 
are in agreement with a straight line through origin of 
slope one. We observe also a deviation from the theoret- 
ical prediction for lower temperatures. This is expected 
since the coherence length approaches the sample length 
for low temperatures, whereas theoretical calculations as- 
sume an infinite system. 

For our experiments, the stability of the effective trap 
depth is crucial, since it defines the reference point of 



4 





1.08- 




1.07 - 




1.06 - 






H 




O 
-to 


1.05- 


o 




II 


1.04- 








1.03 - 








1.02- 




1.01- 




1 - 




50 100 150 
Teal [nK] 



40 60 80 100 120 140 160 180 
Teal [nK] 

FIG. 4: Noise-corrected peak height of the correlation func- 
tion Q2 (x = 0, texp) plotted against the expected temperatures 
corresponding to the calibration T ca i. The vertical error bars 
are given by the statistical fluctuations, the horizontal ones 
indicate the accuracy of the calibration. The theoretical pre- 
diction (grey) is plotted for comparison, its width is given 
by the uncertainty in optical resolution. The inset shows a 
comparison of calibrated temperature T ca i and measured tem- 
perature Tmcas (black squares) to a line through origin with 
unity slope (solid). 



the temperature calibration. The precision of <h(a; 
0. i exp ) is then given by the statistical fluctuation of the 
measurement and the noise characteristics of the imaging 
system. For 300 measurements we achieve an accuracy 
of ~ ± 30 nK, where the main limitation is given by the 
statistical fluctuations. 

To conclude, we have studied the two-point density 
correlation function of quasicondensates in free, collision- 
less expansion. We find an excess of correlations (bunch- 
ing) at short inter-particle distance and an oscillatory be- 
havior at finite distances with correlations below unity, in 
good agreement with the theoretical models of [l3|. We 
show that a quantitative comparison can be used for ther- 
mometry in regimes, where conventional methods based 
on absorption imaging fail. We are convinced that two- 
point density correlations provide a powerful probe for 
low-dimensional systems, nonequilibrium dynamics 
and integrability [19J . 
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For T < 20 nK we would expect a smooth crossover to 
a finite-size condensate; this is however not accessible in 
our setup 

We find it interesting to note, that (in the quasi-lD 
regime) important information on the two-point corre- 
lation function of quasicondensates can be obtained from 
standard absorption imaging without the need for single- 
atom sensitivity or atom shot-noise-limited detection. 
The ability to measure temperatures significantly below 
T c is a unique feature of our setup; it does not invalidate a 
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broad interest in alternative thermometry methods based 
on density correlations. 



